Sloan Digital Sky Survey Quasars in the SWIRE ELAIS Nl Field: 
Properties and Spectral Energy Distributions 
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ABSTRACT 

We present a mid-infrared analysis of 35 quasars with spectroscopic redshifts selected from 
the Spitzer Wide-area InfraRed Extragalactic Survey (SWIRE). We discuss their optical and 
mid-infrared (MIR) colors, and show that these quasars occupy well defined regions in MIR 
color-color space. We examine the issue of type-I AGN candidate selection in detail and propose 
new selection methods based on mid-IR colors. The available multi-band data allows us to 
construct two new, well-sampled quasar templates, covering wavelengths from the ultraviolet to 
the MIR. 

Subject headings: galaxies: active - quasars: general - infrared radiation 



1. Introduction 

Efforts to understand the physics that drive 
the activity in active galactic nuclei (AGN) have 
been underway for decades, and have recently seen 
rapid progress thanks to improvements in both in- 
strument design and theoretical modeling. Obser- 
vationally, new facilities have provided important 
pieces of information, which when combined with 
advances in modeling have allowed a reasonably 
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clear picture of the emission mechanisms in AGN 
to emerge. The optical continuum (up to 1 um) 
of type-I AGN, is seen to be a power law, which 
can be explained by the existence of an accretion 
disk around a central supermassive black hole. At 
longer wavelengths, one starts observing emission 
thought to arise from a dusty torus, which, de- 
pending on the geometry, peaks at wavelengths 
between 15 and 100 pi. This emission is thought 
to be the reprocessed emission of the UV/optical 
radiation from the accretion disk by the particles 
composing the torus, namely silicate and graphite 
grains (e.g. Granato & Danese 1994; Efstathiou 
k Rowan- Robinson 1995; Nenkova et al. 2002). 

Studies with the Infrared Space Observatory 
(ISO) showed that the mid-infrared (MIR) window 
can be used to distinguish AGN from starbursts 
(Rigopoulou et al. 1999; Laurent et al. 2000; Far- 
rah et al. 2002, 2003). Furthermore, recent work 
conducted on a sample consisting of 25 quasars 
detected by the European Large Area ISO Sur- 
vey (ELAIS) at 15 um (Afonso-Luis et al. 2004) 
showed the potential of combined optical and MIR 
observations in the construction of quasar sam- 
ples, and in achieving a better understanding of 
the physical properties of AGN. These studies also 
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showed, however, the need for deeper multi-band 
observations in the MIR in order to refine phys- 
ical constraints on the torus models and unified 
schemes. 

With the advent of Spitzer it is now possible 
to obtain MIR photometry for very large samples 
of galaxies. In particular, the Spitzer Wide-area 
InfraRed Extragalactic Survey (SWIRE; Lonsdale 
ct al. 2003; Lonsdale ct al. 2004) offers an unprece- 
dented opportunity to study the MIR properties 
of AGN thanks to both the large areal coverage 
(~ 50 deg 2 ) and to the number of bands in which 
the fields are being observed (IRAC bands 3.6, 4.5, 
5.8 and 8 pm and MIPS bands 24, 70 and 160 
pm). In this paper we investigate properties of a 
complete optically selected spectroscopically con- 
firmed type-I quasar sample from the Sloan Digital 
Sky Survey (SDSS), consisting of 35 objects. The 
paper is structured as follows. Section 2 describes 
the MIR, optical and near-IR data available for 
this sample; these properties are then used in Sec- 
tion 3 in order to determine type-I AGN candidate 
selection criteria from multi-color data. Section 4 
discusses the SEDs of the objects. Section 5 deals 
with the Black Hole (BH) mass and IR luminosity 
estimates. Notes on individual objects are given in 
Section 6. Finally, Section 7 presents a discussion 
of the results. 

2. The Quasar Sample in the various 
Wavelengths 

The MIR data used here are taken from the 
SWIRE ELAIS Nl field, and were obtained in 
February 2004 with both IRAC and MIPS. For the 
purposes of this work only the four IRAC bands 
(referred to as IRAC1, IRAC2 etc throughout) 
and MIPS 24 um (hereafter MIPS24) are used. 
The SWIRE catalogs we use throughout this work 
were processed by the SWIRE collaboration. De- 
tails about the data can be found in Lonsdale et 
al. (2004), Surace et al. (2004) and Shupc ct al. 
(2004). 

The Sloan Digital Sky Survey has validated 
and made publicly available its Data Release 2 
(DR2; Abazajian et al. 2004), partially covering 
the SWIRE EN1 field. DR2 comprises 32241 
quasars with redshift <2.3 and 3791 quasars with 
redshift >2.3. A total of 35 spectroscopically con- 
firmed quasars lie within the 3.5 deg 2 of SWIRE 



EN1 covered by the SDSS DR2 spectroscopic re- 
lease. Two of the objects only have 24 pm detec- 
tions as they fall in areas at the very border of 
the field and are not covered by IRAC, while four 
others have detections in all four IRAC bands but 
do not have MIPS 24 pm observations due to the 
failure of some MIPS scans during the MIPS cam- 
paign. None of the non-detections are drop-outs, 
instead all are due to the inhomogeneous coverage 
of the field at its edges. Their z-magnitudes reach 
19.1 for objects with redshifts typically less than 
2.3 and go up to a magnitude deeper for higher 
redshifts (Richards et al. 2002). Table 1 shows 
the (optical) positions, redshifts, SDSS AB mag- 
nitudes and MIR fluxes and errors (in mJy) for the 
35 quasars composing our sample. The sequence 
number in the first column in used hereafter as an 
identifier for the objects. 

The quasar sample was matched with the 
2MASS all sky catalog and near-IR counterparts 
were found for six objects, including the highest 
redshift object of the sample at z = 3.653. Six of 
the objects also have ELAIS 15 pm counterparts 
(Vaccari et al. 2004; Rowan-Robinson et al. 2004; 
Afonso-Luis et al. 2004). Details can be found in 
Section 6. 

Fig. 1 (a), (b) and (c) show the optical col- 
ors (in AB) of the 35 quasars (cyan filled cir- 
cles) superposed on the 16710 quasars of the SDSS 
Data Release 1 (DR1; Abazajian et al. 2003), oc- 
cupying the same regions of color space. Fig. 1 
(d) shows the redshift distribution of our sample, 
which peaks at slightly lower redshift as compared 
to the redshift distribution of the DR1 quasar cat- 
alog (Schneider et al. 2003). 

The sample's MIR colors, defined as -2.5 x 
log(flux i /fiux i+ i) with i and i + 1 two consec- 
utive bands, are shown in Fig. 2 (a), (b) and 
(c), superposed on the bandmerged catalog (see 
Surace et al. 2004) extracted from a small region 
(10%, plotted for clarity) of the SWIRE EN1 field. 
The solid, short-dashed, long-dashed and dashed- 
dotted lines correspond to evolutionary tracks of 
a type-I quasar, a red quasar, Arp220 and Sc tem- 
plates respectively (Polletta et al. 2004). The 
quasars' colors are almost independent of their 
redshift due to the essentially power-law nature 
of their SEDs; objects of very different redshift 
can find themselves very close together in this di- 
agram. The Arp220 and Sc templates were chosen 
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because they are those that pass the closest to the 
quasar locus. Quasars are redder than the main 
bulk of the population, at least up to ~ 8 um, 
but then start mixing completely with the rest of 
the objects. Both type-I and red AGN populate 
the redder part of the color diagrams, their very 
red colors being due to torus emission. Fig. 2(d) 
shows the (IRAC1 - IRAC4) color as a function of 
rcdshift and compares it to the evolutionary tracks 
for the various types of objects (linestyle coding 
same as in the other panels of Fig. 2). These 
well defined regions in the MIR color-color space 
can be used standalone criterion for type-I 
AGN candidate selection when optical informa- 
tion is unavailable, and as additional constraints 
to the optical selection and photometric redshift 
determination. 

Models for the IR and submm SEDs of normal, 
star-burst and active galaxies show that quasars 
have, on average, MIR fluxes some 10 to 100 times 
larger than their optical fluxes (Rowan-Robinson 
2001). Fig. 3 compares the r-band flux to the 
IRAC1 flux at 3.6 um (left panel) and MIPS24 flux 
at 24 pm (right panel) for the objects with IRAC 
and MIPS counterparts. The fluxes confirm the 
model predictions, showing that optically bright 
quasars tend to have brighter (M)IR counterparts. 
The solid, dotted and dashed lines denote constant 
IR-to-optical flux ratios, with values 1, 10 and 100, 
respectively. 

3. Quasar Candidate Selection in the MIR 

In this section we discuss type-I AGN candidate 
selection methods. No candidate list is included 
here as the aim is to propose some selection crite- 
ria, describe their advantages and drawbacks and 
the contribution from MIR photometry, and the 
color properties of the candidate sample as a whole 
rather than for individual objects. A similar dis- 
cussion but proposing different selection criteria 
can be found in Lacy et al. (2004), who present 
AGN candidates in the Spitzer First Look Sur- 
vey. For an independent analysis of the location 
of AGN in MIR color space see also Haas et al. 
(2004). 

The MIR band-merged catalog in SWIRE EN1 
was associated with the five-band optical catalog 
derived from the Isaac Newton Telescope Wide 
Field Survey (WFS; McMahon et al. 2001). The 



SWIRE EN1 and INT WFS fields do not have a 
100% overlap; WFS covers only ~ 7 deg 2 of the 
SWIRE EN1 field, down to a limiting AB r-band 
magnitude of 24.1. The final bandmerged catalog 
consists of ~293000 objects, with at least one WFS 
and one Spitzer detection each. Details about the 
optical data and the bandmerged catalogs can be 
found in Surace et al. (2004). 

In optical surveys the use of a morphological 
pre-selection make main sequence stars to be the 
most common contaminants of the quasar candi- 
dates samples. With the addition of MIR data, 
star/quasar separation is straightforward, as stars 
occupy well defined regions in the IRAC color 
space due to the Rayleigh-Jeans regime in their 
SEDs (see Surace et al. 2004 for a detailed dis- 
cussion), separate from the regions populated by 
quasars. In MIR selected quasar samples the ma- 
jor contaminants will instead be normal and star- 
burst galaxies, of both low (z < 0.6) and high (~3) 
redshift (Fig. 2). However, MIR selection can also 
reveal type-II AGN that would mostly be missed 
if only optical selection is used. 

The well defined region our quasar sample oc- 
cupies can be used to constrain the color space 
in which quasar candidates are selected. Fig. 4 
shows the location of the 35 SWIRE-SDSS quasars 
in a (IRAC1 - IRAC2) versus (IRAC2 - IRAC3) 
diagram, superposed with all the sources with an 
r-band counterpart present in the ^7 deg 2 where 
WFS overlaps with the SWIRE EN1 field. The red 
dashed lines delineate the region of quasar candi- 
dates, with equations: 

0.35x(IRACl-IRAC2)+0Alx(IRAC2-IRAC3) = 0.053 

lAx(IRACl-IRAC2)-0.6x(IRAC2-IRAC3) = -0.4 

(IRACA - IRAC2) - (IRAC2 - IRAC'i) = 0.26. 

The region defined here is only indicative, and the 
tolerance depends on the degree of completeness 
and the number of contaminants that are accept- 
able. 

A total of 1870 quasar candidates populate the 
region inside the dashed lines. Ths number drops 
to 1290 when detections in all four IRAC channels, 
MIPS 24 pm and r-band are required. This num- 
ber rises to 2520 when only IRAC channels 1-3 
detections are needed (which are the only relevent 
bands for the selection criterion applied here). 
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signature indicates the presence of nuclear activ- 
ity. 
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Fig. 4.— (IRAC1 - IRAC2) versus (IRAC2 - 
IRAC3) diagram showing the location of the 31 
out of 35 quasars with IRAC1, IRAC2 and IRAC3 
counterparts occupy (cyan filled circles). The 
black contours and points show objects present in 
the ~7 deg 2 where WFS overlaps with the SWIRE 
EN1 field. The red dashed lines delineate the re- 
gion of quasar candidates. For more than 90% 
of all objects the error bars are smaller than that 
shown on the right bottom corner of the figure. 

A type-I AGN sample selected based on MIR 
colors cannot be entirely complete, and could also 
suffer from contamination by spiral and starburst 
galaxies at specific redshifts. When optical mor- 
phology is taken into account the resulting sample 
will be cleaner but less complete, as nearby type-I 
AGN will probably be excluded. Note also that 
the more MIR information one adds the more ac- 
curate the selection can become (see for instance 
Lacy et al. 2004), but then one is limited by the 
substantially lower sensitivity of IRAC channels 3 
and 4. Type-II AGN are more difficult to iden- 
tify. For an in depth discussion on their colors see 
Polletta et al. (2004). Using X-ray selected AGN, 
Franceschini et al. (2004) show that the optical 
to MIR SEDs of type-II AGN can match those 
of normal spiral, starbursting and even passively 
evolving elliptical galaxies, while only their X-ray 



3.1. Model predictions on quasar MIR 
counts 

The number of quasars in this sample, and their 
optical and MIR fluxes can give an estimate of 
the number of AGN SWIRE will detect. All 35 
quasars are at least 3 times brighter than the 
5cr limit at 5.8, 8.0 and 24 pm and at least 25 
times brighter at 3.6 and 4.5 pm, have i < 19.1 
when z < 2.3 and i < 20.1 when z > 2.3. Us- 
ing these limits, expected numbers of AGN have 
been calculated based on a new version of the Xu 
et al. (2003) models, including a new evolution- 
ary model for dusty galaxies. Galaxies with ob- 
scured AGN are those with a MIR excess in the 
rest frame: /as/tm/ houm > 0.2. Their evolution 
function is assumed to be the same as that of op- 
tical quasars, which is a pure luminosity evolution 
function of the form (Boyle et al. 2000): 

L*(z)=L*(0)10 13ez -°- 27z2 (z<7) (1) 

but with L* now defined at 25pm. 

Fig. 5 shows the predicted redshift distribu- 
tions per deg 2 from this model, compared to the 
actual distribution of the 35 SDSS quasars (dashed 
line). The solid line corresponds to the model re- 
alisation using the SWIRE 5tr limits in both IRAC 
and MIPS and the dotted line is the distribution 
predicted when the SDSS optical spectroscopic 
limits are applied as well. Note that the model 
does not distinguish between type-I and type-II 
AGN. The dashed line is the scaled redshift distri- 
bution of the 35 quasars, normalised to 1 deg 2 . 

The total number of objects predicted by the 
model for a 3.5 deg 2 area (which is the coverage of 
SWIRE EN1 by the SDSS DR2) using the SWIRE 
5<7 limits and the SDSS i-band spectroscopic limit 
was 39, i.e. 11 per deg 2 , very close to the actual 
number of quasars in our sample. The predicted 
redshift distribution is also close to the observed 
one, as can been seen by comparing the dotted line 
with the black dashed line in Fig. 5. This good 
agreement is most probably due to the cut in the i- 
band, which most likely excludes all type-II AGN. 
The model predicts some 440 AGN (types I and II 
together) per deg 2 with redshifts up to 4 and de- 
tections in all IRAC bands and the MIPS 24 pm 
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channel down to the SWIRE MIR limits (when 
no optical constraints are imposed), with some 20 
z > 3 objects per deg 2 . The difference between 
the observed and estimated numbers implies a 
large amount of incompleteness introduced in the 
quasar samples by optical constraints. The num- 
ber of objects (1870 in the ~7 square degrees of 
SWIRE EN1 covered by the WFS data) previously 
selected as quasar candidates is 1.6 times smaller 
than the one predicted by the model. This number 
drops to 1.2 (i.e. 360 candidates per deg 2 ) when 
no optical counterpart is required. The model es- 
timates for the AGN differential number counts 
(N/deg 2 /dex) at 3.6 and 24 pm for various optical 
magnitude limits are shown in Fig. 6. Note that 
on these plots only the respective MIR detections 
are required, increasing the estimated number by 
a factor of almost 3.5 with respect to the previous 
model realisation where all IRAC band and MIPS 
24 pm are required. With very few exceptions, 
the quasars in our sample have much brighter 3.6 
and 24 pm fluxes than the fluxes where the model 
counts peak. 

The missing population consists of type-I 
quasars omitted from the optical selection, type- 
I AGN without all IRAC and MIPS detections 
as well as type-II AGN that are mostly optically 
obscured or optically extended. In fact, Lacy et 
al. (2004) suggest that up to 50% of their MIR 
selected AGN could be sufficiently obscured not 
to be observed by SDSS. 

4. Spectral Energy Distributions 



Fig. 5. — Histograms showing the number of AGN 
per deg 2 expected for the SWIRE 5a limits (solid 
line), the number of AGN expected when the 
SDSS spectroscopic limits are applied (dotted line) 
and the actual quasars (dashed line). 



The Spitzer data in SWIRE EN1 fill a large 
wavelength gap in the SEDs of AGN. The quasar 
sample discussed here is optically selected and 
complete down to i = 19.1. Therefore, when con- 
sidering their SEDs, one does not see the behav- 
ior of an IR selected sample but rather the IR 
properties of optically selected objects with MIR 
counterparts. Figs. 7 and 8 show the observed 
SEDs (with 2MASS and ELAIS 15 pm measure- 
ments included whenever available; see Section 6 
for details) in AF\ (erg/sec/cm 2 ) overplottcd on a 
template constructed from optically selected (from 
the Palomar-Green (PG) Survey; hereafter quasar 
template) type-I AGN (from Polletta et al. 2004). 
Observations and templates have been normalised 
at 3.6 pm (or to the r-band in the few cases where 
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the 3.6 pm flux was not available). The quasar 
template was built by extending the composite 
SDSS quasar optical spectrum (Vanden Berk et 
al. 2001) into the IR with the average SED of 14 
PG quasars. At A < 0.7 pm the template is the 
SDSS composite quasar spectrum; the averaged 
PG template is used only at A > 0.7 um. 

The vast majority of the objects have a concave 
SED, with a change of slope at around 1 um in the 
rcstframe. The few exceptions will be discussed in 
detail in Section 6. The general shape of the SEDs 
roughly follows the quasar template but the use of 
all 35 quasars provides better sampling, and allows 
for the refinement of the template, especially in the 
IR. 

Using the 35 SWIRE-SDSS quasars, two new 
templates were constructed. The first one (dashed 
line in Fig. 9) is the average of all of the data, nor- 
malised at 0.3 um in 18 wavelength bins between 
0.1 and 23 um (diamonds). The upper template 
(dotted line) is the average of the data points taken 
from the 25% objects with the highest infrared 
emission in the same wavelength bins (triangles). 
The number of flux values in each bin varies from 
7 to 37 with an average of 20 when all sources are 
used, and from 3 to 10 with an average of six in 
the case of the second template. In the optical, at 
A < 0.6 pm for the first template and A < 0.3 pm 
for the second, the more detailed quasar template 
was used since it fits well the average data of both 
templates. At A > 20 pm, since no data are avail- 
able, we use the quasar template after scaling it to 
match the two new templates. The solid curve is 
the original quasar template shown in Figs. 7 and 
8, shown here for comparison. A detailed study of 
the individual SEDs, tori model template fitting 
and parameter analysis is outside the scope of this 
work and will be presented in a separate paper. 
The two new IR quasar templates defined above 
are shown in Tables 2 and 3. 1 

Fig. 10 presents the bolometric luminosities 
(computed assuming an (fi, A, H ) = (0.3, 0.7, 70) 
cosmology), integrated from ^800 A to 24 pm, as 
a function of redshift, revealing 85% of them to 
be ultraluminous (with L^ > 10 12 L ) and one 
hyperluminous {Lboi > 10 13 L Q ). The bolomet- 
ric luminosities have been computed assuming the 



The SEDs can be downloaded also from the url 
http://www.iac.es/proyect/swire 



newly derived SEDs. The Scott effect clearly dom- 
inates this flux-limited sample as more distant ob- 
jects need to be intrinsically brighter in order to 
be part of it. 

5. Black Hole Masses 

In this section we investigate a possible correla- 
tion between bolometric luminosity and BH mass, 
as the relationship between these two quantities 
is a measure of the Eddington ratio. Under the 
assumption that the dynamics of the Broad Line 
Region (BLR) are dominated by the gravity of 
the central supermassive BH, the latter can be es- 
timated from the relation Mbh — RblrV 2 /G, 
where Rblr is the radius of the BLR and V is 
the velocity of the line-emitting gas. Tradition- 
ally, V is estimated from the full width at half 
maximum (FWHM) of the H^ line in emission 
(Kaspi et al. 2000). However, for quasars with 
redshifts greater than z ~ 0.8, is no longer 
seen in observed-frame optical spectra, and Mg 
II has been suggested as an alternative estimator 
(McLure & Jarvis 2002). 

Here the BH mass is computed as suggested by 
McLure & Dunlop (2004): 

M»L = 2 f A£ 3ooo \ 0,62 ( FWHM(MgII) \ 2 
M Q ' \10 37 W ) V km/ sec ) 

(2) 

and 

Mbh / AL 5100 \ 061 ( FWHM{H p ) \ 2 

M e \10 3 ?Wj \ km/ sec J (> 

Since the redshift range in which Mgll can 
be observed with optical telescopes is large we 
chose to use the relation between Mgll and the 
(monochromatic) continuum luminosity at 3000 
A (Eq. 2) whenever available and only use H3 
(Eq. 3) in the few cases when Mgll is not present. 
We make use of the values of the emission lines 
as measured by the SDSS pipeline and limit the 
study to objects with redshifts lower than, typi- 
cally, z = 2.1 as at higher redshift both estima- 
tors fall outside the observed spectral range. Fig. 
11(a) shows the BH mass distribution as a func- 
tion of redshift, with objects for which Mgll (H^) 
was used for the estimation plotted in filled circles 
(open squares). The error bars are estimated from 
the uncertainty of the FWHM of the lines. Black 



6 



hole masses tend to increase with redshift up to 
z ~ 1 but then they stabilise around a value of 

5 xlO 8 M . The very low BH mass value of the 
object at z ~ 2 comes from the apparently un- 
derestimated value of the Mgll line, due to a very 
noisy spectrum. An attempt to measure this line 
manually gave a value three times larger for the 
Mg II line, increasing the value of the BH mass by 
roughly an order of magnitude (open symbol in 
Fig. 11). For the rest of the objects, the manual 
measurements of the emission lines agree within 
2a with the ones coming from the SDSS pipeline. 
For the six objects lying in the redshift range from 
0.5 to 0.83 for which both Mgll and EL3 were avail- 
able, BH masses were computed using both esti- 
mators. For four of the objects the difference in 
the masses were less than a factor of 1.3, for the 
other two, however, the estimated masses were al- 
most one order of magnitude apart, with H^ giving 
the lowest value in both cases. 

Fig. 11(b) shows the evolution of the bolomet- 
ric luminosity as a function of BH mass. A trend 
is seen up to a bolometric luminosity of ~ 10 12 
L Q but it is of low significance, with a Spearman 
correlation coefficient of 0.32. This is due to small 
number statistics. Note that the determination of 
BH mass, dependent as it is on optical luminosity, 
can only be applied to unobscured sources. 

6. Notes on individual Objects 

Four of the quasars have radio detections: 
source 18 at 1.4 and 4.85 GHz (~ 80 and 171 
mJy, respectively; Gregory & Condon 1991; White 

6 Becker 1992); and sources 4, 19 and 27 with 
FIRST 20cm detections and fluxes 4.20 ± 0.15 
166.23 ± 0.146 mJy and 43.82 ± 0.148 mJy, re- 
spectively. Furthermore, two of the radio sources, 
18 and 19, also have X-ray counterparts from the 
ROSAT All Sky Survey. Finally, six objects have 
ISO 15 urn fluxes and six have 2MASS J, H, K 
counterparts, as shown in Table 4, along with the 
radion information. 

Object 18, at a redshift of 0.245, is particu- 
larly interesting. As can be seen in Fig. 7 (sec- 
ond object in the last row), its SED does not have 
the concave shape that the majority of the objects 
have. Its SED is instead convex, and can be fit- 
ted by two power laws, a relatively flat one in the 
optical and a much steeper one in the IR. Even 



though this object is a type-I quasar with broad 
emission lines, as can be seen in the left panel of 
Fig. 12, the behavior of its SED is that of a blazar. 
Its inverted radio spectrum (higher flux in higher 
frequencies) also points towards this conclusion. 

A similar behavior is observed for source 22 at 
z = 0.214, whose spectrum is shown in the right 
column of Fig. 12. Unfortunately no other data 
are available for this object. Both objects show 
signs of variability, an important characteristic of 
all blazars, between the times the spectroscopy 
and the photometry were made, as can be seen in 
the lower panels of Fig. 12, where the spectra and 
magnitudes of the two objects are shown in the 
AB system. For a comprehensive study of quasar 
variability based on this type of comparison see 
Vanden Berk et al. (2004). The photometric er- 
rors are of the size of the filled circles (or smaller) 
and can not be seen on the plot. 

Finally, sources 21 and 24, whose optical spec- 
tra are considerably redder than the quasar tem- 
plate (Fig. 8), are in fact broad absorption line 
quasars, as visual inspection of its spectrum re- 
vealed. Recent studies relate reddened quasars 
with the presence of broad absorpion lines in their 
spectra (e.g. Reichard ct al. 2003 and references 
therein). 

7. Discussion 

This paper presents an analysis of a sample 
of 35 spectroscopically confirmed SDSS quasars 
within the SWIRE EN1 field, observed with 
Spitzer as part of the SWIRE Legacy Survey. 
Even though SWIRE EN1 is a small part of the 
whole 50 deg 2 SWIRE survey (-18%) and SDSS 
DR2 only covers a third of it, this study gives a 
good initial indication of the AGN science that 
can be carried out with SWIRE data. 

All SDSS quasars within the SWIRE EN1 field 
have Spitzer counterparts. The distinct location 
of the spectroscopically confirmed quasars on MIR 
color-color plots greatly contributes to the selec- 
tion of typc-I AGN candidates, independent of 
their optical properties. They are redder than the 
majority of the population at least up to 8 pm 
(restframe from 1.7 to 6.1 pm for the redshift range 
covered) and continue forming a rather compact 
clump in the redder colors. Appropriate color- 
color combinations can be found in order to sepa- 
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rate quasars of specific (but large) redshift ranges 
from the various types of contaminant galaxies. 
A comparison between model predictions and se- 
lected quasar candidates indicates the existance of 
a large number of obscured AGN, some of which 
probably reside among the proposed candidates. 

BH masses, estimated from emission lines, show 
a tendency with redshift up to z ~ 1, however 
the small number statistics do not allow any con- 
clusion on the dependence or not on bolomctric 
luminosity. 

Even though large multi-wavelength spectro- 
scopic efforts are still required in order to un- 
derstand the nature and propeties of AGN, we 
demonstrate that Spitzer MIR data contribute sig- 
nificantly towards this goal. 
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the Jet Propulsion Laboratory, California Insti- 
tute of Technology under NASA contract 1407. 
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Fig. 1. — Optical (AB) colors (panels a, b and c) and redshift distribution (panel d) of the quasar sample 
(cyan filled circles). Black contours and black points indicate the location of the ~ 16700 spectroscopically 
confirmed quasars of SDSS DR1. 
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Fig. 2. — MIR colors of the quasar sample (filled cyan circles). The black points are the objects in the 
bandmerged catalog of tile 3_1 (see text for details). The black solid, red short-dashed, green long-dashed 
and blue dashed-dotted evolutionary tracks correspond to a type-I quasar, a red quasar, Arp220 and Sc 
templates, respectively (from Polletta et al. 2004). For more than 90% of all objects the error bars are 
smaller than that shown on the right bottom corner of panel (a). 
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Fig. 3. — 3.6 pm (left panel) and MIPS 24 pm (right panel) versus optical (r-band) flux shown for the 
objects with IRAC1 and MIPS 24 pm counterparts, respectively. The solid, dotted and dashed lines denote 
constant IR-to-optical flux ratios, with values 1, 10 and 100, respectively. The errors on the IR fluxes are 
smaller than the plotting symbols and therefore cannot be seen on the figures. 
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Fig. 6. — AGN number counts for various optical (r-band) limits predicted, at 3.6 (left) and 24 pm (right) 
fluxes. Only the respective MIR detections are required, increasing the estimated number by a factor of 
almost 3.5 with respect to Fig. 5. 
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Fig. 7. — Observed SEDs for the 35 SWIRE-SDSS quasars (rilled circles) overplotted to an average quasar 
template from a subsample of the PG survey described in the text, shifted to the redshift of the object. The 
objects appear in the same order as in Table 1, i.e. ordered by RA. 
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Fig. 8.— Observed SED for the SWIRE-SDSS quasars (continued). 
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Fig. 9. — Average (dashed line) and higher MIR flux (dotted line) templates derived from the SWIRE-SDSS 
quasars. Diamonds and triangles show averaged data points in various wavelength bins. The solid line is 
a template derived from the PG sample making use of the SDSS composite spectrum in the optical and is 
showed here for comparison. The grey circles are the observed data points. 
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Fig. 10. — Bolometric luminosities, integrated from ^800 A to 24 pm, as a function of redshift expressed in 
units of solar luminosity, computed assuming the new SEDs. 
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Fig. 11. — (a) BH masses versus redshift. In open squares (filled circles) objects for which (Mgll) was 
used for the estimation, (b) BH mass as a function of bolometric luminosity. The open circle indicates the 
value of the BH mass of the dubious object as measured by hand. 
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Fig. 12. — Upper panels: SDSS spectrum of objects 18 and 22 showing broad emission lines, characteristic 
of a typc-1 quasar. Object 22 also presents strong narrow forbidden lines, as in type-II AGN. Lower panels: 
spectra and magnitudes of the two objects are shown in the AB system, demonstrating signs of variability. 



18 



Table 1 

Positions, redshift, optical magnitudes and aperture fluxes for the 35 spectroscopically 
confirmed quasars. the sdss magnitudes are in the ab system, fluxes are given in mjy. 
Note that the errors in the IRAC1 and IRAC2 fluxes are lower than 0.01 mJy (typically 

OF THE ORDER OF A FEW JiJy) AND ARE NOT REPORTED HERE. 
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.50 


1 


.65 


2, 


.09 


± 





.01 


2 


.87 


± 


0, 


.01 


13.8 


± 


0.02 


23 d 


16: 


:07 


:05. 


,16 


53: 


:35: 


:58 


.5 


3 


.653 


23 


.61 


19 


.35 


18. 


.15 


18 


.14 


17 


.98 





.40 





.40 





.66 


± 





.01 


1 


.33 


± 


0. 


.01 


5.89 


± 


0.02 


24 


16: 


:08 


:28. 


.32 


53 


:52: 


:51 


.9 


2 


.960 


25. 


,83 


21 


.25 


20. 


.44 


19 


.80 


19 


.62 





.12 





.11 


0. 


.14 


± 





.01 





.17 


± 


0, 


.01 


0.61 


± 


0.02 


25 


16: 


:08 


:56. 


,78 


54 


:03: 


:13 


.7 


1 


.917 


19, 


.29 


19 


.25 


19, 


.23 


19 


.28 


18 


.94 





.20 





.25 


0. 


.39 


± 





.01 





.60 


± 


0. 


.01 


1.34 


± 


0.02 


26° 


16: 


:09 


:08. 


.95 


53 


:31: 


:53 


.2 





.816 


19, 


,05 


18 


.50 


18. 


.48 


18 


.68 


18 


.49 





.66 





.78 


1. 


.13 


± 





.01 


1 


.44 


± 


0, 


,01 








27 e 


16: 


:09: 


:13. 


,18 


53 


:54: 


:29 


.5 





.992 


18. 


,63 


18 


.34 


18. 


,04 


18 


.00 


18 


.01 





.62 





.68 


0, 


.96 


± 





.01 


1 


.34 


± 


0. 


,01 


5.65 


± 


0.03 


28° 


16: 


:09: 


:43. 


,67 


53 


:30: 


:41. 


.0 


1 


.328 


19. 


,42 


19 


.01 


18 


,54 


18 


.39 


18 


.30 





.59 





.98 


1. 


.37 


± 





.01 


2 


.00 


± 


0, 


.01 








29° 


16: 


:09: 


:50. 


,72 


53: 


:29: 


:09 


.5 


1 


.716 


18. 


.15 


18 


.04 


18. 


,04 


17 


.86 


17 


.79 





.58 





.97 


1. 


.74 


± 





.01 


2 


.89 


± 


0. 


.01 








30 a < f 


16: 


:09: 


:53. 


,02 


53: 


:34: 


:43 


.8 


1 


.212 


18. 


.10 


17 


.96 


17. 


,83 


17 


.89 


17 


.97 





.69 


1 


.06 


1. 


.75 


± 





.01 


2 


.49 


± 


0. 


.01 








31 c 


16: 


:10: 


:07. 


.11 


53: 


:58: 


:14 


.0 


2 


.030 


18. 


.99 


18 


.92 


18, 


.85 


18 


.77 


18 


.56 





.19 





.26 


0. 


,52 


± 





.01 





.93 


± 


0. 


.01 


3.26 


± 


0.03 



Table 1 — Continued 



Seq Nr RA Dec z u g r i z IRAC1 IRAC2 IRAC3 IRAC4 MIPS 24 

32 16:10:57.70 53:19:46.2 1.576 19.61 19.42 19.25 18.99 19.07 0.17 0.25 0.48 ± 0.01 0.78 ± 0.01 2.35 ± 0.03 

33 16:12:15.12 53:12:47.7 1.300 19.61 19.39 18.95 18.89 18.86 0.59 0.98 1.45 ± 0.01 2.16 ± 0.01 4.38 ± 0.02 
34 d 16:12:38.27 53:22:55.0 2.138 17.92 17.82 17.82 17.72 17.48 0.41 0.49 0.82 ± 0.01 1.30 ± 0.01 3.56 ± 0.02 
35 16:13:38.23 53:13:51.5 1.304 18.49 18.53 18.27 18.29 18.41 0.32 0.45 0.68 ± 0.01 1.03 ± 0.01 3.72 ± 0.02 



a Object in the MIPS 24 pm gap 

b Object at the border of the field not (or partially) covered by IRAC 

c Object with ELAIS 15 pm detection 

d Object with 2MASS JHK counterparts 

e Object with radio detection 

f Object with IRAS detections 

s Object with RASS detections 

c,d ' e 2MASS magnitudes, ISO 15 pm and radio fluxes can be found in Table 4 



Table 2 

Quasar SED template created using the 35 SWIRE-SDSS quasars, spanning the 

WAVELENGTH RANGE BETWEEN ~ 800 AND 2 X 10 6 A. THE FLUXES ARE NORMALISED AT A = 0.3 
THE COMPLETE VERSION OF THIS TABLE IS IN THE ELECTRONIC EDITION OF THE JOURNAL. T] 

PRINTED EDITION CONTAINS ONLY A SAMPLE. 



A (jim) fA/fo.3A» A (P m ) V f o.3 M A (pm) f A /fo.3 M A (P m ) V f o.3 M A (pm) f A /fo.3 M 



0. 


.079 


0, 


.369 





.195 


1 


.904 





.311 





,966 


0, 


,428 





.489 


0.544 





.293 


0. 


083 


0, 


.546 





.199 


1 


.778 





.315 


0. 


.979 


0. 


,432 





,541 


0.548 





.294 


0. 


087 


1. 


.346 





.203 


1 


.691 





.319 


0. 


.969 


0. 


,436 





.539 


0.552 





.298 


0. 


091 


2. 


197 





.207 


1 


.741 





.323 


0. 


,944 


0. 


,440 





.481 


0.556 





.309 


0. 


095 


3, 


190 





.211 


1 


.611 





.327 


0, 


,908 


0. 


,444 





.455 


0.560 





.313 


0. 


099 


3, 


,293 





.215 


1 


.513 





.331 


0, 


,876 





,448 





.457 


0.564 





.304 


0. 


103 


4, 


.694 





.219 


1 


.446 


o 


.335 


0. 


.852 


0. 


.452 


o 


.466 


0.568 





.290 


0. 


.107 


3 


.550 





.223 


1 


.413 





.339 


0. 


.853 


0. 


,456 





.465 


0.572 





.278 


0. 


111 


3 


128 





.227 


1 


.439 





.343 


0. 


,840 


0. 


,460 





.450 


0.576 





.274 


0, 


115 


3, 


.006 





.231 


1 


.495 





.347 


0. 


,807 





,464 





,440 


0.580 





.284 


0, 


119 


4, 


,272 





.235 


1 


.444 





.351 


0, 


,771 


0. 


,468 





,417 


0.584 





.317 


0. 


.123 


7, 


,052 





.239 


1 


.436 





.355 


0, 


,743 


0. 


,472 





.398 


0.588 





.382 


0. 


127 


3 


,706 





.243 


1 


.476 





.359 


0. 


,729 





,476 





.391 


0.843 





.210 


0. 


131 


3, 


,405 





.247 


1 


.432 





.363 


0. 


,711 





,480 





,426 


0.995 





.159 


0. 


135 


3, 


,038 





.251 


1 


.381 





.367 


0, 


,693 


0. 


,484 





,534 


1.505 





111 


0. 


139 


3 


,644 





.255 


1 


.324 





.372 


0, 


690 


0. 


,488 





.534 


2.276 





.094 


0, 


143 


2. 


,804 





.259 


1 


.270 





.376 


0. 


,665 


0. 


,492 





.455 


3.442 





.078 


0. 


.147 


2, 


,667 





.263 


1 


.238 





.380 


0. 


,624 


0. 


,496 





.449 


5.205 





.049 


0. 


.151 


2, 


,910 





.267 


1 


.197 





.384 





,620 


0. 


,500 





.467 


7.871 





.025 


0. 


.155 


4, 


,357 





.271 


1 


.227 





.388 


0, 


,607 


0. 


,504 





.431 


11.90 





.017 


0. 


.159 


2, 


.695 





.275 


1 


.315 





.392 


0. 


,571 


0. 


,508 





.339 


18.15 





.014 


0. 


163 


2 


,667 





.279 


1 


.874 





.396 


0, 


,552 


0, 


,512 





.339 


23.14 





.010 


0. 


.167 


2. 


,427 





.283 


1 


.284 





.400 


0. 


,533 


0, 


,516 





.352 


29.51 





.007 


0. 


171 


2. 


,296 





.287 


1 


.141 





,404 


0. 


,514 


0. 


,520 





.351 


37.63 





.005 


0. 


175 


2. 


,245 





.291 


1 


.067 





,408 


0. 


,531 


0. 


,524 





.348 


47.98 





.003 


0. 


179 


2. 


,126 





.295 


1 


.087 





.412 


0. 


,523 


0. 


,528 





.346 


61.17 





.002 


0, 


183 


2. 


,059 





.300 


1 


.000 





,416 


0. 


,490 


0. 


,532 





.336 


77.99 





.001 


0. 


.187 


2. 


,235 





.303 





.940 





,420 


0, 


,477 


0. 


,536 





.315 








0. 


191 


2. 


,562 





.307 





.954 





.424 


0. 


.472 


0. 


.540 





.300 
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Table 3 

Quasar SED template created using the SWIRE-SDSS quasar subsample with the largest 
IR fluxes. The complete version of this table is in the electronic edition of the Journal. 

The printed edition contains only a sample. 



A (11m) fA/f .3 M A (ym) f\/f .3» A (ym) A (ym) ix/fo.^ 






.079 





.369 


0. 


148 


2 


.674 


0. 


,216 


1 


.514 


0.285 


1 


,204 


0. 


.082 





.645 


0. 


150 


2 


.829 


0. 


,219 


1 


.446 


0.288 


1 


,132 





.084 





.850 


0. 


153 


4 


.000 


0. 


,222 


1 


.436 


0.291 


1 


,060 





.087 


1 


.418 


0. 


156 


3 


.251 





.225 


1 


.419 


0.294 


1. 


,098 





.090 


2 


.027 


0. 


159 


2 


.695 


0. 


,228 


1 


.448 


0.296 


1 


,076 


0. 


.093 


2 


.975 


0. 


162 


2 


.663 


0. 


.231 


1 


.494 


0.300 


1, 


.000 


0. 


096 


2 


.931 


0. 


165 


2 


.540 


0. 


.233 


1 


.487 


0.417 


0. 


.601 





.099 


3 


.293 


0. 


168 


2 


.368 





.236 


1 


.433 


0.513 


0. 


.526 


0. 


.102 


3 


.579 


0. 


170 


2 


.289 





,239 


1 


.436 


0.631 


0. 


.487 





.105 


3 


.455 


0. 


173 


2 


.275 


0. 


,242 


1 


.505 


0.776 





.380 


0. 


.107 


3 


.488 


0. 


176 


2 


.185 


0. 


,245 


1 


.430 


0.955 


0. 


.265 


0. 


.110 


3 


.073 


0. 


179 


2 


.126 


0. 


,248 


1 


.429 


1.037 


0. 


.234 





.113 


3 


.018 


0. 


182 


2 


.059 


0. 


,251 


1 


.395 


1.929 


0. 


.180 





.116 


3 


.123 


0. 


185 


2 


.199 


0. 


,253 


1 


.360 


3.587 


0. 


,117 





.119 


4 


.272 


0. 


188 


2 


.294 





.256 


1 


.307 


6.671 


0. 


.064 


0. 


.122 


9 


.076 


0. 


190 


2 


.717 


0. 


,259 


1 


.270 


12.40 


0. 


,031 


0. 


.125 


5 


.074 


0. 


193 


2 


.035 


0. 


.262 


1 


.256 


21.00 





.024 


0. 


.127 


3 


.421 


0. 


196 


1 


.821 


0. 


.265 


1 


.173 


30.23 


0, 


,014 


0. 


.130 


3 


.525 


0. 


199 


1 


.778 





,268 


1 


.185 


43.53 





.008 


0. 


.133 


3 


.282 


0. 


202 


1 


.722 





,271 


1 


.216 


62.67 


0, 


,004 


0. 


.136 


3 


.041 


0. 


205 


1 


.711 


0. 


,274 


1 


.279 


90.24 





.002 


0. 


.139 


3 


.644 


0. 


208 


1 


.735 





.276 


1 


.345 


129.9 


0, 


,001 


0. 


.142 


2 


.940 


0. 


211 


1 


.644 


0. 


,279 


1 


.874 








0. 


.145 


2 


.702 


0. 


213 


1 


.559 


0. 


,282 


1 


.347 
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Table 4 

2MASS Vega magnitudes, ISO fluxes (from Vaccari et al. 2004) radio fluxes and related 

ERRORS (WHENEVER AVAILABLE) FOR THE QUASARS WITH NEAR-IR AND/OR 15 pM AND/OR RADIO 

COUNTERPARTS. FOR REFERENCES SEE TEXT. 



Seq Nr 




RA 






Dec 








./ 






H 






K 






ISO 15 pm 
(mJy) 


4.85 GHz 
(mJy) 


1.4 GHz 
(mJy) 


A 


1 Pv 


.oy .oo. 


1 3 

. 10 


04±, 


■Z19-IT? 


.o 


10. 


01 


_i_ 

nz 


n 


90 


i a 7Q _i_ 




1 A 


Q9 -4- 


n 


10 






a on 4- n i k 

4:.ZU in U.±0 


1 

Iz 


10 


■no. kc\ 


,yo 


04, 


OU.O 1 . 


n 

y 
























i.yy ± u.oo 






17 


16: 


:05:23. 


.10 


54: 


:56:13. 


.3 
























2.02 ± 0.39 






18 


16: 


:05:46. 


11 


54: 


:39:09. 


.8 


16. 


.90 


± 


0. 


.18 


16.18 ± 0. 


.22 


15 


.57 ± 


0. 


.18 




171 


80 


19 


16: 


:06:23. 


.56 


54: 


:05:55. 


.8 


16. 


.30 


± 


0. 


.14 


16.14 ± 




15. 


.37 ± 


0. 


.22 


3.91 ± 0.63 




166.23 ± 0.146 


20 


16 


:06:30. 


.60 


54: 


:20:07. 


.5 
























3.14 ± 0.53 






21 


16: 


:06:37. 


.88 


53: 


:50:08. 


.4 
























1.46 ± 0.28 






22 


16: 


:06:55. 


.34 


53: 


:40:16 


.8 


16. 


.37 


± 


0. 


10 


15.33 ± 0. 


10 


14. 


.32 ± 


0. 


.07 








23 


16: 


:07:05. 


.16 


53: 


:35:58. 


.5 


16. 


.60 


± 


0. 


.12 


16.29 ± 0. 


22 


15 


.72 ± 


0. 


.22 








27 


16: 


:09:13. 


.18 


53: 


:54:29. 


.5 




























43.82 ± 0.148 


31 


16: 


:10:07. 


.11 


53: 


:58:14. 


.0 
























1.59 ± 0.30 






34 


16: 


:12:38. 


.27 


53: 


:22:55. 


.0 


16. 


.70 


± 


0. 


15 


16.15 ± 0. 


15 


15 


.54 ± 


0. 


23 









